Generalized momenta of mass and their 
applications to the flow of compressible fluid 



Olga Rozanova 

Mech.&Math. Faculty, Moscow State University, Moscow, 119992, Russia 
Abstract 

We present a technique that allows to obtain certain results in the compressible 
fluid theory: in particular, it is a nonexistence result for the highly decreasing at 
infinity solutions to the Navier-Stokes equations, the construction of the solutions 
with uniform deformation and the study of behavior of the boundary of a material 
volume of liquid. 
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1 Main equations 



The motion of compressible viscous, lieat-conductive, polytropic fluid in M x 
M", n > 1, is governed by the compressible Navier-Stokes (NS) equations 

9tp + div,(pt;) = 0, (1.1) 

dt{pv) + div^(pf ®v) + V^p = DivT, (1.2) 

dt (^pI^P + Pe) + div^. (i^p\v\^ + pe + p)v^ = div(Tt;) + kAJ, (1.3) 

where p, u = {ui, ...,Un), p, e, 9 denote the density, velocity, pressure, internal 
energy and absolute temperature, respectively; T is the stress tensor given 
by the Newton law T = Tij = p{diVj + djVi) + Xdivu6ij, with constant 
coefficients of viscosity p and A (/i > 0, A + -/i > 0), A; > is the coefficient of 
heat conduction. We denote Div and div the divergency of tensor and vector. 
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respectively. The state equations are p = (7 — l)pe, and p — Rp9, where 
7 > 1 and R > are the specific heat ratio and the universal gas constant, 
respectively. Thus, we can consider (NS) as a system for unknown p, u, p. 
Indeed, (NS) and the state equations give 

n k V 

dtp + {v, V^p) + 7p divv = (7 - 1) E '^iA^i + -5 (1-4) 

For classical solutions (NS) is equivalent to system (1.1, 1.2, 1.4), denoted 
(NS*) for short. 

(NS*) is supplemented with the initial data 

li p — X — k — 0, we get the gas dynamic (GD) equations. 



2 Conservation laws, generalized momenta of mass and the decay 
rate 



System (NS) is the differential form of conservation laws for the material vol- 
ume V{t) (i.e. the volume, consisting from the same particles). (NS) expresses 
conservation of mass m — J pdx, momentum P — f pvdx, and total en- 

V(t) V(t) 

ergy S = J (^p\v\'^ + pe) dx = Ek{t) + Ei{t). Here Ek{t) and Ei[t) are the 

v(t) ^ ^ 
kinetic and internal components of energy, respectively. 

When V{t) — R", the conservation of mass, angular momentum and energy 
take place provided the components of solution decrease at infinity sufficiently 
quickly. 

Definition 2.1 We say that the classical solution solution to (NS) belongs to 
the class K{M{t), a), if there exist a positive vector- function M{t) — {M^{t), Mp{t), 
Mp(t)) with components from C([0, 00)), a constant vector a = {a^, ajj,^, a^. ap, ae) 
and constants Ro > 0, T > 0, such that \v{t,x{t))\ < M^(t)|x(t)|"% \Dv{t,x{t))\ < 
MB„(t)|a:(t)|"°", p{t,x{t)) < Mp(^)|a;(^)|"^ p{t,x{t)) < Mp{t)\x{t)\"^ , 9{t,x{t)) < 
Me(t)|a;(t)|"« for all trajectories x{t) such that \x{t)\ > Rq, t > T. 

Let us introduce a functional 

G4>{t) = J p{t,x)(j){\x\)dx. 

V(t) 
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If = then G^{t) is the usual momentum of mass. By analogy for 

others (f) we call G^{t) the generalized momentum of mass. 



2.1 Decay rate for (NS) equations 



Let us choose the appropriate class K[M{t).,a) to guarantee conservation of 
mass, momentum, energy and convergence of the mass momentum G{t) for 

(NS) system. It is sufficient to set a = ctp, «p, q;^) = {—n,—n — 

1, —n — 2 — e, —n — £, —n) with a constant £ > 0. We denote this class K^s- 

If the heat conductivity is zero, we do not need to require the decay of 9, i.e 
a = {—n, — n — 1, —n — 2 — e, —n — £, ae) with an arbitrary last component. 
We denote this class K^So- 



2.2 Decay rate for (GD) equations 



In this case the behavior of velocity is less restrictive. Here it is sufficient to 
set a — (q;^,, aovi ~ 2 — £, —n — e, ae) with e > 0, < 1 and arbitrary 
a^v and ae. The components of velocity may rise as |x| ^ oo. We denote this 
class Kgd- 



3 Behavior of generalized momenta of mass on solutions 



We denote a = {ai, ...,(7k) the vector with components = ViXj — VjXi, 
i> j, ij = l,...,n, k = 1,...,K, K = C^. 

The following lemma describes the behavior of generalized momenta for (GD) 
system. 

Lemma 3.1 Let us suppose that (l){\x\) belongs to the class inside V{t). 
For classical solution to system (GD) the following equalities take place: 

G',{t)= f ^^{v,x)pdx, 

V(t) 



where 



J \X\ J \X\ 

V[t) V{t) 
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V(t) dV{t) 

where v is a unit outer normal to V. 

The proof is a direct application of the Stokes formula. 

Let us consider 0(|a;|) = \\x\'^ and V{t) = M". We denote the respective 
functional as G{t). In this case the following lemma is true. 

Lemma 3.2 For classical solutions to (NS), k = 0, of the class K^So ^.i^'d 
solutions to (GD) of the class Kcd the following equality holds: 

G"{t) = 2Ek{t) + n{^-l)E,{t). (3.1) 

PROOF. First of all, in (NS) case G"{t) = h{t) + hit) + h{t) + h{t) + h{t) , 
where Ii{t) = J ^^^^pdx, hit) = J ^-^pdx, h{t) = n J pdx, hA^) = 

— lim / {x,h')pddBji, hit) = J(J)ivT,x)dx = lim / TijXiVj ddB^ — 

{2fi + nX) J divvdx = hm / {TijXi — {2fi + nX)vi6ij)h'j ddBji, where -Bj? : 

\x\ < R, we sum the repeated indices. Owing to respective assumptions 
on the decay of solution as |x| — >■ cx) all improper integrals converge and 
h{t) = hit) = 0. To end the proof we note that hit) + h{t) = '2Ek{t), h{t) = 
n{'~f — l)Ei{t). In (GD) case the integral h is missing. 

Remark 3.1 Equality (3.1) for (GD) was obtained firstly in Jl^]. 

Remark 3.2 It follows from (3.1) that n{-f -1)8 < G"{t) < 28. The inequal- 
ity implies the same two sided estimate of G{t) that was obtained in 3D for 
the sum of squared semi-axis of a non-rotating gas ellipsoid next to vacuum 
If. 



4 Necessary conditions for existence of global solution 



If the initial density po is compact, then in arbitrary space dimensions no 

+2, exists [3|. This blowup 



solution to (NS) from G\[0, oo), //"(R")), m > 
result depends crucially on the assumption about compactness of support 
of the initial density. Thus, the question remains: is it true that the global 
in t smooth solution exists for any smooth initial data in the case where 
the support of initial density coincides with the whole space? Below we find 
necessary conditions of existence of the global solution with prescribed decay 



rate as x 



oo. 
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Theorem 4.1 If the global int solution to (NS) of the class K^Sq exists, then 
the solution components grow as t — > oo at least such that 



t 

j M^{t) dr = 0{t^-"-), Mp{t) = 0{t^+'), e > 0, t ^ oo. (4.1) 





PROOF. It follows from Lemma 2 that G"{t) > (7 - l)£, 

Git) > llziM^s ^ ^ ^(Q^_ ^4 2) 

Let us get the estimate of Git) from above. 

We consider a material volume V(t) that initially coincides with the ball |x| < 
i?o (see Definition 1). We denote -Bi?(t) = {x ■.\x\< R{t)) a ball that contains 
V{t). Prom the definition of K^s we have 



dt 

Further, 



< R'^^it) M,{t), R{t) < (1 - a,) J M^(t) dr + R^ 



1— a„ 



G{t)^l i I p\x\'dx+ J 

W(t) R»\B,i 

^Mp{t) J \x\'^+'^'' dx < 



ax | < — - — m+ 




M,{t) 1 1 M,{T)dT 



If the growth of Mi,(t) and Afp(t) is less then prescribed in (4.1), then the 
latter inequality contradicts to (4.2). 

Remcirk 4.1 In particulary, Theorem 1 implies that there exist no smooth 
solution to (NS) from the class K^So ^^^^ My{t) — const and/or Mp{t) — 
const. 

The analog of Theorem 1 in the case of (GD) equations is the following. 

Theorem 4.2 // the global in t solution to (GD) from the class Kgd ex- 
ists then the functions Mp{t) are M^it), restricting the density and velocity 
behavior as t ^ 00, grow at least such that Mpit) = Oit"^'^"), e > 0, and 
/q M„(r) dr = O(lnt) (for — 1) or as prescribed in (4-6) (for <1). 
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The statement can be proved exactly as Theorem 1. 



5 Motion with uniform deformation 

We dwell in this section on (GD) equations, however due to the specific choice 
of the velocity field the results remain true for (NS) with k = 0. 

The motion with uniform deformation (i.e with linear profile of velocity v{t, x) = 
A{t) X, A{t) is a matrix (n x n)) was considered in many works. Let us men- 
tion [1], where this special solution to (GD) was firstly constructed in 3D and 
[S], where many applications are given. Generally speaking, for this solution 
the velocity, density and pressure may blow up as t — T < oo and |x| — > oo. 
Below we show how the solution with uniform deformation can be constructed 
by means of the generalized momenta of mass. Moreover, the requirement of 
finiteness of energy and mass momentum prohibits the blow up, therefore the 
respective solution is globally in time smooth. 

Below we compare the procedure of constructing of the simplest solution with 
the linear velocity profile 

v{t, x) = a{t) X (5.1) 
for a different choice of generalized momenta. 

5. 1 Basing on the momentum of mass 

Let us consider once more the particular case 0(|a;|) = the respective 

functional is the usual momentum of mass, G{t). 

Lemma 2 and (1.4) imply 

G'{t) = 2a{t)G{t), a'{t) = -a\t) + KG-'''^^ {t), (5.2) 

with a constant > depending on initial data. If the initial data satisfy the 
compatibility condition 

Po(kl) = -(7 - l)G-i(0)i?,(0)po(|x|), (5.3) 

then the density and pressure can be found from (1.1) and (1.4) as 

t t 
p(t, |x|) = exp(— n j a(T)dT)pQ{\x\exp{— J a(T)dT)), (5.4) 
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t t 
|x|) = exp(— n7 j a{T)dT)po{\x\exp{~ J a{T)dT)). (5.5) 



It follows from (5.2) that a{t) solves 

a\t) = -a^(t) + Ki exp ^^-((7 - l)n + 2) J a(r) rfrj , (5.6) 

with Ki = K(G(0))-«T-i)"+2)/2. Since Ki > 0, the solution a{t) remains 
bounded for all t > and a(t) = 0(t~^), t — > 00 [f\. 

Remark 5.1 In the case of globally smooth solution discussed above, the ve- 
locity v(t,x) = a(t)x, therefore = 1. Here M^it) = a{t) = 0{t~^), and 
Jq M^{t) dr = O(lnt), Mp{t) = 0(t^+''), t ^ 00. Thus, according Theorem 1, 
the solution presents a less possible rate of its component in t. 

Remark 5.2 In the solutions with linear profile of velocity for (GD) equa- 
tions (including the presence of the Coriolis force and damping) were con- 
structed for a general matrix A{t) in 2D and for particular cases in 3D. 



5.2 The "excluding pressure" case 



Let us choose as (/) a function, proportional to the fundamental solution of the 
Laplace operator, namely 

0(|x|) = ln|a;|, n = 2, 0(|x|) = n ^ 2. 

In this case Lemma 1 implies 

G'^{t) = X,{n) l^^pdx, 

R" ' ' R" ' ' ^ ' 

where 

Ai(?T-) = 1, 72 = 1,2 and 2 — n, n > 3, 

X2{n) = 0, n = 1; —1, n = 2, and {l — n){2 — n), n>3, 

X^ln) =0, n = 1; 1, n = 2, and 2 — n, n > 3. 

Here we take as V(t) the space M", all improper integrals are supposed conver- 
gent both at the origin and at infinity. Actually, this signifies that v{t, 0) = 0. 

We use the value of pressure only in the origin; in this sense the pressure in 
the remaining space is excluded. 
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Let us consider the case n > 3. Here G'^(t) = \i{n) a{t)G^{t), 

Glit) = X^in) a'{t) G^it) + ^^"'"j exp (-^n f a(r) dr 
The latter system imphes 

a'(t) = -a'(t) + /r2exp I -((7- l)n + 2) / a(r) I , (5.7) 



-a'(t) + /r2exp ^-((7-l)n + 2) ^ a(r) drj 



P(0 0)G"~'^ (0) 

With = We can see that (5.6) coincides with (5.7), the only 

difference is in constants Ki and K2. Provided initial data satisfy the compat- 
ibility condition 

GMv'om = 4^^Po(|x|)|x|, (5.8) 

the density and pressure can be found by formulas (5.4), (5.5). 

The analysis of equation (5.6) (or (5.7)) shows that if the constant Ki (or 
K2) is positive, as in our case, then the solution exists for all t > and 
a{t) = 0(t~^) as t 00. 

We see that in the second case (0 = |xp~") we get a more wide class of 
solutions. Indeed, in Section 5.1 we need to require the significant decay of 
solutions as |x| — CX3 to guarantee the finiteness both of energy and momentum 
of mass, i.e. we restrict ourselves by solutions from Kgd with ay = 1. When 
we construct solutions using the generalized momentum, we do not require 
the finiteness of energy and the behavior of smooth pressure and density, 
prescribed by convergence of G^{t) and condition (5.8) is the following: 

p(|x|) = 0(|xn, g<0, p{\x\) = 0{\x\'), s = q-2<0, \x\ ^ 00. 



6 Behavior of boundary of a liquid volume 



We mention here very briefly one more application of the generalized momenta 
of mass to (GD) equations (see [S] for details). Namely, the expansion of 
boundary of a material volume inside of a smooth flow of gas can be studied 
by this method. The last question is connected with a problem of air pollution: 
one can be interested if under usual meteorological conditions a polluted cloud 
will attain some geographical object. 

We suppose that initially a point xq do not belong to V(t). Let us set the 
following question: what conditions we have to impose on initial data provided 
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they are known only inside of V(0), to guarantee that the boundary of given 
material volume within a smooth flow will attain a given e - neighborhood of 
point Xq? It is clear that for the answer to this question certain assumptions on 
the thermodynamic variables in the whole space have to be done. To formulate 
these assumptions we need the following definition. 

Definition 6.1 We say that the pressure in the moment t is distributed along 
the boundary dV (t) of domainV{t), xo ^ V(t), regularly with constant M > 0, 



dV{t) ^' ^ 



< M, where u is a unit outer normal. 



If p{t,x) is constant, then the integral in the left hand side of the latter in- 
equality is equal to zero. Setting M sufficiently small we assume that the 
material volume will not meet a zone of large gradient of pressure. 

To prove the following theorem we use the generalized momentum of mass 
with = \x — Xq\'^, q < 0. 

Theorem 6.1 [8] Let a finite material volume V{t) of compressible liquid with 
- smooth boundary dV{t) do not contain a point xq. Suppose that the flow 
is - smooth for all t E [0, T], T < oo, and the pressure along the boundary 
dV{t) is distributed regularly with a constant M uniformly in t for t G [0,T]. 

Let us choose some real numbers q < — and e, < £ < dist(9V(0), xq). 
Then for all initial data there exists such constant 5 < 0, depending on initial 
data, £, g, T, M, n, 7 that if initially 

X — Xq\'^~'^ (^(0, x),x — Xq) Po{x) dx < 5, 



V(0) 

then within a time ti, later then T, the boundary of given material volume will 
attain the e - neighborhood of point xq. 
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